arXiv: 1501.07903vl [physics.optics] 30 Jan 2015 
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We show experimental and numerical evidence of spontaneous self-symmetrization of focused laser 
beams experiencing multi-filament at ion in air. The symmetrization effect is observed as the multiple 
filaments generated prior to focus approach the focal volume. This phenomenon is attributed to the 
nonlinear interactions amongst the different parts of the beam mediated by the optical Kerr effect, 
which leads to a symmetric redistribution of the wave vectors even when the beam consists of a 
bundle of many filaments. 


I. INTRODUCTION 

An intense laser pulse propagating in air undergoes 
various self induced transformations if its peak power ex¬ 
ceeds a critical value [T] Per ~ 4.7 GW. The optical Kerr 
effect is responsible for an increase of the refractive index 
leading to beam self focusing, a cumulative process evolv¬ 
ing towards a collapse. An electron plasma is efficiently 
generated via optical field ionization when the beam in¬ 
tensity has increased above ^ 10^^ W/cm^. Ionization is 
associated with nonlinear absorption of laser energy and 
with a decrease of the local refractive index. Both ef¬ 
fects act against the growth of intensity by self-focusing 
and eventually arrest the collapse. The interplay between 
nonlinear effects leads to the formation of a narrow light 
filament leaving a plasma channel in its wake, surrounded 
by a laser energy reservoir. This reservoir maintains an 
energy flux toward the filament core that compensates 
for nonlinear absorption. This optical entity can live over 
several meters and be generated at long distances 
Ionization in the atmosphere by filamentation has been 
reported at a distance of 1 km from the laser [4]. When 
the reservoir is exhausted and fails to feed the filament 
efficiently, a slowly diverging bright channel is observed 
and can be detected at several kilometers away from the 
laser source. Several teams have reported an effect called 
beam self-cleaning in the low power regime P ~ Per [S]- 
[9]. Moll and Gaeta showed that slightly elliptical beams 
with power close to Per tend to increase their roundness 
and recover a circularly symmetric profile when they un¬ 
dergo self-focusing, just before ionization. This effect has 
been interpreted as due to beam collapse and reshaping 
into a universal self-similar spatial profile, the Townes 
mode um, which is intrinsically spatially symmetric. 
Prade et al. [6] measured conical emission for ultraviolet 
pulses undergoing filamentation in air and reported self- 
improvement of the beams’s spatial mode quality, sug¬ 
gesting the generation of propagation invariant modes 
called nonlinear X-waves [TTl |T2]. Liu and Chin have 
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shown that the beam cleaning process only occurs for 
the frequency downshifted components of the supercon¬ 
tinuum generated during the filamentation process. They 
attributed this effect to the role of self-focusing as a spa¬ 
tial filter, filtering out the high-order mode components 
of the beam whereas the fundamental mode is focused, 
producing a high quality filament core. Naturally this 
process is promoted for redshifted frequencies which are 
generated in the leading edge of the pulse and undergo 
mainly self-focusing and diffraction [8]. From these ex¬ 
amples, it is clear that even if beam self-cleaning received 
different interpretations, this effect occurs for beam pow¬ 
ers moderately above the critical power. 

For P ^ Per, however, a beam rises multiple filaments 
across the profile, as manifestations of the dynamics fol¬ 
lowing local collapse. These filaments are born on beam 
intensity fluctuations [nmu and grow due to modula¬ 
tion instability [151 |T6] (^I)- Several filament interactions 
manifested as cross phase modulation (XPM), fusion, re¬ 
pulsion, spiraling, and fission have been reported in this 
high power regime [isiiTTHin]. It has also been shown 
recently that focused multi-filamented beams maintain 
a substructure during the focal region and may produce 
thick filaments associated to unusually high intensity and 
plasma levels, coined as superfilaments [21]. Up to date, 
all studies done with high power beams focus on the fila¬ 
ment features, rather than on the whole beam dynamics. 

In this work we investigate an effect of the global 
structure of powerful (P ^ Per) beam profiles that re¬ 
mained unexplored. We demonstrate experimentally and 
by means of numerical simulations that under focusing 
conditions, powerful beams undergo self-induced sym¬ 
metrization in the global scale above certain threshold 
power, Pfh ^ Per, even if hundreds of filaments are born 
along propagation. As a consequence, the beam washes 
out all interactions and distortions it might have suffered 
prior to the focus. However the nature of this very robust 
phenomenon differs from the previously reported beam 
self-cleaning effect associated with a single filament, in 
the sense that a fundamental mode does not necessar¬ 
ily emerge after self-induced symmetrization. We used 
two indicators to distinguish between beam self-cleaning 
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and self-induced symmetrization: the global beam sym¬ 
metry degree ^ and the beam quality factor . For the 
single filamentation regime, beam self-cleaning reported 
previously is essentially associated with an improvement 
of the beam quality factor, even if the beam symmetry 
degree may improve as well. In contrast, in the case of 
high power beams considered here, only global symmetry 
improves substantially, whereas the beam quality factor 
always remains far from that of a Gaussian beam or the 
Townes mode (M^ ^1). To illustrate this effect, we 
performed experiments and numerical simulations on the 
nonlinear propagation of heavily distorted input beams 
characterized by a dramatic decrease of their symmetry 
degree after passing through intensity masks. Results 
clearly demonstrate that the combination of the Kerr ef¬ 
fect and focusing is crucial to observe self-symmetrization 
of highly asymmetric beams. 

Self-symmetrization of powerful (multiterawatt) laser 
beams may prove useful when it is important to homo¬ 
geneously illuminate a target placed in the focal region. 
This effect also allows us to obtain a relatively uniform 
cylindrical plasma channel around focus which shape is 
independent from the quality of the initial laser beam 
profile. Such property is particularly interesting for ap¬ 
plications based on laser filamentation at high power such 
as guiding of electric discharges [22ll24] and contact-less 
capture of currents [25], control of aerodynamic flows 
[2^ . or lasing effect in air [22MI]. 


II. EXPERIMENTAL AND NUMERICAL 
PROCEDURES 


A. Experiments 


The collimated output 50 fs pulse from a linearly po¬ 
larized multi-Terawatt Ti:Sa CPA laser (Enstamobile) at 
Ao = 800 nm is weakly focused with an / = 5 m lens in 
air at atmospheric pressure. Burn patterns from the laser 
beam are recorded on calibrated presolarized Kodak pho¬ 
tographic plates at different distances from the focusing 
lens. This simple technique reliably detects in a single 
shot the stage of evolution of multiple filaments formed 
out of a single laser beam [2]. Plasma channels give rise 
to characteristic ^ 50—100 /im wide circular burns on the 
photographic plate, which are easily distinguished from 
the ^ 1 mm circular dark spots from bright channels 
and the less intense energy reservoir. Similar measure¬ 
ments were performed by intentionally adding chirp to 
the pulse to adjust the peak intensity. The same self- 
symmetrization effect of multiple filaments was observed 
for the transform limited 50 fs pulses and for pre-chirped 
pulses (see figure captions). This allows us to dismiss 
temporal effects as the origin of symmetrization. 


B. Beam characterization 

For each recorded beam intensity profile, we computed 
the center of mass (CM) on the XY- (transverse) plane 
and then considered the radial intensity traces (with fixed 
polar angle 0), /^(r) (r = 0 at the CM). We define the 
degree of symmetry, able to resolve the internal struc¬ 
ture (intensity fluctuations) of the beam profiles, as the 
weighted average of the intensity fluctuations: 0 < ^ < 1 
(see appendix) 


r)-I<pAr)} 
r) + I<pAr)} 

( 1 ) 

which measures the similarity in between all pairs of the 
radial intensity traces, {/^(r),/^/(r)}. Perfect circularly 
symmetric profiles are characterized by ^ = 1, whereas 
beams with a large asymmetry are characterized by a 
relatively low symmetry degree. 

Regarding the beam quality factor, M^, we use the 
widely used definition by Potemkin et al. [32| for a beam 
with complex field f (x, y) = A(x, y) exp(i0) and intensity 
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and r = (x^y). Since the calculation of the quality factor 
requires knowledge not only of the beam intensity but 
also of the spatial phase, we characterized beams in terms 


of only for the results of our numerical simulations. 



C. Simulations 

Propagation of the monochromatic field envelope 
f (x,^, z) in air (no ~ 1), with frequency ujq = 27rc/Ao, is 
modeled by means of a unidirectional beam propagation 
equation accounting for diffraction, optical Kerr effect, 
multiphoton absorption, plasma absorption, and plasma 
defocusing, respectively [33] : 

^ _ 

dz 

^ £+i‘^n2\£\'^8-\ (/38|^r^ + <J[1 + iujoTc]p) £. 
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FIG. 1. Experimental transverse profiles of an undistorted 
focused beam. Images recorded on the photographic papers 
at various distances for input pulses of 50 fs (not pre-chirped), 
carrying 130 mJ each and having a beam width of ~ 3.5 cm. 
The lens has its geometrical focus at 5 m. Beam and pulse 
widths are given at the intensity FWHM and propagation 
distances are marked on top of each image. The 1 cm length 
bar is included for reference. 


Here n 2 = 2 x cmVW, /^g = 8 x 10“^^ 

cr = 5.6 X 10“^^ cm^, and Tc ~ 350 fs. Our (2 + 1)D 
modeling does not contain temporal dynamics. Hence, 
plasma effects are accounted for in the so-called frozen 
time approximation where the explicit standard rate 
equations describing electron generation by optical field 
ionization and avalanche are used to generate a mapping 
of the electric field peak intensity X = \£\‘^ to the electron- 
plasma density, p{X). This calculation is performed for 
referenee pre-chirped Gaussian pulses with peak inten¬ 
sity X with duration given by experimental conditions. 
The electron density, p(X), used in the propagation model 
Eq. 0 is determined from this mapping, at the tempo¬ 
ral center of the reference pulse. This procedure provides, 
as shown previously and below, numerical results in good 
agreement with experiments uni HD EH- Simulations are 
initialized with 10% intensity and 0.2% phase noise in 
order to mimic experimental irregularities in the input 
beam. 


III. RESULTS 
A. Undistorted beams 

In the absence of any mask, a beam with power P ^ 
550 Per undergoes typical multi-filamentation dynamics. 
Figure shows the darkening pattern impressed by such 
a beam at several distances from the focusing lens on 
photographic plates. Already after 2 m of propagation 
several ionizing filaments are present. These converg¬ 
ing filaments grow in number and connect into networks 
upon further propagation (> 4 m). Before the geometric 
focus, the ionized zone fills a central circular spot of a 
few mm in diameter and ^ 1 m long, where filaments 
are in close contact ED. The number of ionized spots is 
significantly reduced to a few which are located near the 
center of the beam. Losses incurred by going through 
the focus have been measured with a calorimeter placed 
before and after the focus. 82% of the initial beam en¬ 
ergy is found in the beam at 3.5 m after the focus. This 
rather high value could be connected to the idea that in 
certain situations the filaments are connected stronger to 
plasma defocusing than to plasma absorption [IT]- Be¬ 
low we show that independently on the intensity mask 
used to distort the input beam profile, the beam recov¬ 
ers after focus the same symmetry degree, hence similar 
roundness, as the one shown in the final stage of Fig. 1. 
In contrast, low energy beams (P ^ Per) do not exhibit 
any improvement of their symmetry degree along prop¬ 
agation since the final beam profile corresponds to the 
diffraction pattern induced by the mask, with identical 
symmetry degree as the initial distorted beam. 

B. Self induced symmetrization of highly distorted 
beams 

In order to visualize the self symmetrization effect we 
have put different masks on the path of the input beam, 
just before the focusing lens (at 0 m). Examples are 
shown in Figs. and where a quarter circle (or Pac- 
man), square, half-plane, and slit masks are used. Sim¬ 
ilarly to the undistorted case, 89% of the input chirped 
pulse energy was measured 3 m beyond the focus. In¬ 
spection of the burnt papers in Figs. and a) reveals 
that ionizing filaments get first organized along patterns 
dictated by the profile of the input beam m- Indeed, 
close to the sharp intensity jumps created by the mask, 
large intensity modulations drive filamentation in a de¬ 
terministic manner. Upon further beam propagation, the 
self-symmetrization effect appears around the geometric 
focus symmetrizing the slowly diffracting output beam, 
despite its initial strong distortion. Solid experimental 
evidence of the symmetrization effect is provided by com¬ 
paring the profiles in Fig. 2 recorded at 1.5 m before and 
after the focus (i.e., at 3.5 and 6.5 m). 

A detailed comparison between experimental and nu¬ 
merical results is presented in Fig. [^for the case of the 
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FIG. 2. Experimental beam profile evolution under spatial 
distortion of the pre-chirped 700 fs pulses carrying 295 mJ 
before the masks (P/Pcr ~ 10^). The different masks screen 
(a) 25 %, (b) 10 %, and (c) 50 % of the input energy. Input 
beam widths and focusing conditions are those of Fig. 
Propagation distances are indicated on top of each column of 
images. 
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FIG. 3. Beam evolution after distorting the input beam with 
a 1 cm wide slit mask, screening ~ 64 % of the input energy 
~ 300 mJ. (a) Experimental and (b) simulated beam inten¬ 
sities along propagation at selected distances, marked in (a). 
Pulse width is 100 fs (pre-chirped) and energy ~ 100 mJ: 
P/Pcr ~ 210 (after mask). Simulations are initialized with 
a peak intensity of ~ 250 GW/cm^. Position of geometrical 
focus is at 5 m and intensities given by the color bar are in 
Bels. 


highest input distortion. Here, a 10 mm wide slit trans¬ 
forms the input circular profile into a rectangle with as¬ 
pect ratio ^ 1/4. Upon propagation the symmetriza- 
tion degree increases as the beam approaches and goes 
through the focal region at ^ 5 m. Figure |^a) shows 
symmetrization of the beam along propagation for the 
beam profiles presented in Figs. [^-b. For a comparison, 
two additional numerical results for linear propagation 
and for an input intensity /q = 5 GW/cm^ are also pre¬ 
sented. The symmetry degree is seen to be higher by the 
end of the nonlinear propagations, reaching its maximum 
after the focus and holding it for a distance of at least 
^ 2 — 3 m. The minimum of the symmetry is observed 
for the linear propagation at the focal plane. This is sim¬ 
ply due to the sharp horizontal line {y = 0) apparent in 
the beam profiles at this position, 2 : = 5 m (see below, 
e.g., in Fig. 6), corresponding to the Fourier transform of 
the slit induced sharp profiles along x (which is a sinc(x) 
function at focus). One can easily see that the width of 
the Fourier transforms of a Gaussian, (at FWHM), 
and a slit, (full width of the cetral ”sinc” lobe), sat¬ 
isfy Aks/Akg = w/Ls X 7r/\/ln2, where w ^ 3.5 cm is 
the FWHM of the input {z = 0) Gaussian profile along 
y and Lg = 1 cm the slit width (along x). It therefore 
follows that whilst Lg ^ w, the horizontal line appear¬ 
ing at focus will be much longer than the spot size along 
y affecting substantially the symmetry degree (see the 
sharp minimum ^ « 0.5 in Fig. 4a). This feature disap¬ 
pears gradually as the input peak intensity is increased 
in simulations up to Jq ~ 5 GW/cm^, because of the 
novel spatial frequencies that are generated (see discus¬ 
sion in section V). Also, the photographic paper used in 
experiments has a different definition in the low intensity 
parts of the beam and therefore the threshold intensity 
at which the horizontal line disappears is expected to be 
different than for simulations. Note in the linear prop¬ 
agations for the masked beams (not shown) the beam 
profiles at a fixed distance before and after the focus are 
almost equal to each other after a rotation of tt in the 
XY plane. The equality does not hold exactly due to the 
input noise and mask induced fluctuations, however the 
difference is small at sight. This is the reason why ^(z) 
is (almost) symmetric from the focal plane (see Fig. 4a). 
In our experiments, the slit mask is the one screening a 
biggest fraction (^ 64%) of the input 300 mJ carried by 
laser pulses, which was the maximum available energy. 
We believe this is why the experimental symmetry de¬ 
gree along propagation (circles in Fig. 4a) starts to drop 
a bit 2 — 3 m after focus. This feature may be appreciated 
even visually in Fig. 3(a), conversely to what is shown in 
Fig. 2 for the other masks. Additionally, the symmetry 
degree presents a local minimum around focus, which is 
a reminiscence of the linear propagation (indeed we see 
similar ^{z) trends in our monochromatic modeling for 
^0 ^ 10 GW/cm^). Still, the symmetry degree is remark¬ 
ably high during 2-3 m and we expect (see below) that 
this would improve for higher input pulse energies. 

The quantitative agreement in between our monochro- 
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FIG. 4. (a,b) Symmetry degree of the beam profiles versus propagation and input peak intensity, respectively, (c) Smoothness 
along propagation for Townes mode (dashed), undistorted focused Gaussian beam propagating in the linear regime (solid-red), 
and heavily distorted beams, corresponding to the linear and Iq — 250 GW/cm^ cases shown in (a). Results shown here 
correspond to the slit mask (see also figures 3 and 5). 


matic numerical simulations and experimental results in 
Figs. l^andQa) strongly suggests that the spontaneous 
symmetrization shown here is an effect essentially dom¬ 
inated by spatial dynamics of the beam. We have com¬ 
puted the symmetry degree at a fixed distance z = 9 m 
(4 m after focus) and varying input power (Fig. 4(b)). 
These measurements reveal an increase of ^{P/Pcr) ex¬ 
hibiting saturation: ^{P/Pcr ^ 40) ^ 0.95. Such be¬ 
havior vividly manifests the nonlinear (intensity depen¬ 
dent) nature of the spontaneous symmetrization. Indeed, 
a systematic numerical study reveals that the Kerr ef¬ 
fect play s a major role in self-symmetrization (see sec¬ 
tion below). We also characterized the beam quality 
via beam quality factor [32l [35]. A Gaussian beam at 
waist is characterized by = 1 and is considered as a 
reference for perfect quality. Figure 4(c) shows that the 
highly distorted beams always present a beam quality 3 
to 5 orders of magnitude worse than that of a Gaussian 
beam at focus or the Townes mode. Therefore, even if 
there is a partial improvement of with propagation 
distance, the main effect is self-symmetrization. Beam 
self-cleaning occurs in the sense of an improvement of 
the beam quality factor but the latter effect is not so ef¬ 
ficient as for beam self-cleaning obtained in the context 
of lower power beams [siiaiH], where PjPcr ^ 1- 


IV. NUMERICAL ANALYSIS OF KEY 
PHYSICAL EFFECTS IN BEAM 
SYMMETRIZATION 

A. Impact of the different nonlinear effects 

We consider below the case of an input beam with a 
peak intensity of Jq = 50 GW/cm^ under the distortion 
of the slit mask. This situation corresponds to the mini¬ 
mum peak intensity needed to observe a high symmetry 
degree (see Fig. 4(b)) and only a few filaments may be 
observed on the output beam profile, note there are only 
four filaments in Fig. 5 (marked by the squares). In order 
to show that the main responsible for self-symmetrization 


of powerful focused beams is indeed the optical Kerr ef¬ 
fect, we made a set of simulations with identical initial 
conditions as those in Fig. 5 in which the different non¬ 
linear terms in Eq. (1) are switched on or off at will (see 
Fig. 6). Below, we are only showing beam profiles at and 
after the focus {z > 5) because those before the focus are 
indistinguishable by sight from those in Fig. 3 (b) for all 
cases. 


When the only nonlinear effect is MPA (n 2 , cr = 0, 
Fig. [^a)), propagation does not differ substantially from 
the one observed in the linear regime, except for the big 
absorption induced across the profile that lead to the ap¬ 
pearance of vertical fringes. Linear propagation would 
lead to the observation of the diffraction pattern of the 
slit. Nonlinear absorption plays the role of a distributed 
stopper, leading to modulations in the diffraction pat¬ 
tern in an effect similar to the Arago spot effect. Nonlin¬ 
ear absorption participates to the fringe formation since 
it enhances the self-healing process of a beam that is 
known to reshape Gaussian beams into Bessel-like beams 
j364[38] . Addition of plasma absorption and defocusing 
(cr ^ 0, Fig. I^b)) has a strong impact in the profile 
after focus: generation of plasma leads to refractive in¬ 
dex jumps along the sharp intensity edges and the de¬ 
focusing induces strong scattering towards the directions 
perpendicular to those edges (light propagates towards 
the higher index regions). As a consequence of this, the 
scattering of light after focus occurs predominantly at an 
angle 7r/2 from the long axis of the slit. It is only when 
the Kerr effect is switched on that symmetrization degree 
improves substantially, as shown in Fig. ic) (/38 = 0) 
and Fig. Note from Fig. id) that with Kerr effect 
symmetry is higher when both MPA and plasma effects 
are accounted for, presumably due to the fact that MPA 
and plasma tend to induce scattering along perpendicular 
directions (at least in the case of the slit mask). Remark¬ 
ably, the symmetrization effects persists even for beams 
experiencing a large multifilamentation (see Fig. 3). 
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FIG. 5. Beam symmtrization in the low filament limit. Loga¬ 
rithmic scale plots of beam Intensity profiles at various prop¬ 
agation distances (see labels) for a beam focused with a lens 
of 5 m of focal length and perturbed with a 10 x 35 mm^ slit 
mask (as in Figs. 3 and 4). Input peak intensity /o = 50 
GW/cm^. All nonlinear effects of Eq. 1 are switched on (n 2 , 
/^s, cr 7^ 0). 


B. Impact of diffraction: non focused beams 

The second essential ingredient for symmetrization (to¬ 
gether with Kerr) is the focusing induced by the lens. We 
have done extensive numerical modeling with collimated 
beams and controlled the density of filaments via the in¬ 
put power, P. None of these simulations showed traces 
of spontaneous symmetrization in the distances of ^ 10 
m. Figure shows an example with the slit mask. Here, 
in order to decrease filament separation, the input power 
was 4 times larger than for the focused cases shown in 
Fig. 3(b) or Fig. 4 with the largest /. The multiple 
filaments develop along propagation and the phase gra¬ 
dients impinged by the mask make some of the optical 
energy migrate into the initially unpopulated area. How¬ 
ever, symmetry is seen to remain rather low along the 
10 m, which implies that the self-symmetrizing effect re¬ 
ported here is intimately related to the focusing induced 
interaction amongst the different parts of the beam. 
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FIG. 6. Influence of the different nonlinear terms, (a)-(c) 
Beam prohles at z = 5,8 m from numerical results corre¬ 
sponding to simulations with identical conditions as in Fig. 

(a) 772, cr = 0, (b) 772 = 0, (c) CT = 0. Golor scale is that 
of Fig. (d) Evolution of symmetry degree for (a)-(c) and 
Fig. HI The nonlinear effects that are on in each simulation 
are listed in the legend. 


V. DISCUSSION: A POSSIBLE EXPLANATION 
FOR THE SYMMETRIZATION EFFECT 

In this section we attempt to give a plausible explana¬ 
tion for the symmetrization effect on the basis on the re¬ 
sults reported and described above. First of all we notice 
that results in Sec. IV show that symmetrization occurs 
mainly due to the Kerr term and that it is not necessar¬ 
ily linked to the presence of many filaments in the beam 
profile. This strongly suggests that symmetrization oc¬ 
curs as a consequence of the Kerr-induced generation and 
annihilation of spatial frequencies (this is indeed a neces¬ 
sary condition). Along this line it is well known (see e.g., 
[39] . for the description in time domain) that the Kerr 
term induces generation of spatial frequencies more effi¬ 


ciently in the regions where intensity gradients are bigger 
(e.g., sharp edges impinged by the mask on the beam). 
The transverse wavenumber of these waves satisfies the 
proportionality relation oc d±I {d± denotes spatial 
derivative in some direction on the transverse plane XY). 
Note here two things: (i) the new spectral component is 
generated at a sharp edge and travels inwards towards 
the intense part of the beam, and (ii) this wave being 
generated, implies that other frequencies present at that 
edge have been annihilated, i.e., the Kerr effect boosts 
energy from the strongly divergent spatial frequencies at 
an edge to new frequencies traveling in the opposite direc¬ 
tion over the transverse plane (this is of course equivalent 
to the qualitative argument explaining why the Towns 
mode [To] exists). In the case of the slit mask, the image 
described above plays a central role along x (horizontal) 
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powerful beams. The absence of symmetrization for col¬ 
limated beams highlights the importance of focusing in 
the global organization in the form of circularly symmet¬ 
ric beam. A possible explanation for this effect is also 
provided. 
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APPENDIX: SYMMETRY DEGREE 


FIG. 7. Simulated beam intensity cross sections along prop¬ 
agation for a collimated cw beam under distortion of the slit 
mask. Absence of beam symmetrization and of filament den¬ 
sity increase are evident (compare with Fig. |^b)). Numbers 
in figures mark propagated distances along z. The reference 
Gaussian pulse used in this modeling is of 100 fs and the in¬ 
put beam peak intensity is 2.5 TW/cm^. All figures have a 
cross-section of 4 x 4 cm^. 


and not so much along y (vertical) because the Gaussian 
profile is much smoother than the rectangular shape. As 
a consequence of this, one could say that Kerr acts as 
a moderator which kills waves with high \k±\ and there¬ 
fore symmetrizes the spectral distribution. However, this 
is not enough to explain the effect reported in this pa¬ 
per. We believe that the second necessary condition for 
symmetrization is the presence of the external focusing. 
Not only because it enhances intensity levels and the ef¬ 
ficiency of the Kerr effect, but also because it tends to 
localize all the frequency components of the beam into 
a reduced focal region that acts as a quasi-point source 
with symmetric spectrum in the transverse plane. This 
unavoidably leads to an after-focus beam that exhibits a 
substantially high overall roundness and we can therefore 
say that the beam self-symmetrized. A precise answer to 
why this effect still happens in the multi filamentation 
regime (P > 150 Per) and what is the role of the tur¬ 
bulent, by nature, multi filament dynamics m requires 
further work and understanding. 


First of all, for a given 2D beam profile, we locate the 
center of intensity (center of mass) 


{x,y)cM 


+ yuy)l{x,y) 


(8) 


at which we locate the origin of the polar coordinates, 
(r,0): r = 0 at {x^y)cM- Then we consider the radial 
intensity traces at a fixed polar angle, /^(r), and then we 
compare them by pairs. Given a pair of these traces we 
define the local relative difference 


D{r, (1)1, (p 2 ) 


UAr) - hAr) 

hAO + UAO 


UAr) - U^jr) r_. 


(9) 

Since D values are more significant in regions where the 
average intensity, (70^,02 (^))5 locally high, we define 
the weighted average relative difference in between the 
pair of traces: 


WQ{(j)i,(f>2) 


IrXD 


W 


e [0,1], 


( 10 ) 


which nullifies for equal intensity traces and tends to 1 
for very dissimilar ones. The overall symmetry degree 
is obtained by comparing all possible pairs of intensity 
traces: 


VI. CONCLUSION 

We have demonstrated a self-symmetrization effect in 
air occurring for high power laser beams experiencing 
multifilamentation and the action of an external focus¬ 
ing force. The collective organization process is very 
different from self-cleaning and is linked to the Kerr 
induced isotropic wavenumber redistribution leading to 
beams with improved overall symmetry. Simulations are 
in excellent qualitative and quantitative agreement with 
the observed features of this effect induced by focusing 


2 r27r n(f)i 

$ = i-^y^ #1 d(f>2WD{(f>i,(f>2)e[o,i], (11) 

where the meaning of ^ is reversed from that of i.e., 
^ is good (bad) for values close to 1 (0) and the nor¬ 
malization factor is d(j)i (i02 = 27r^. Eqns. (9-11) 
trivially combine to give the symmetry degree presented 
in Eq. (1). In the numerical analysis we have used the 
discrete version of the above equation taking into account 
only a finite number of intensity traces ( 0 i ,2 ^ m,n), 
which adopts the form: 












2 M-1 M 

^ = ^ ^ «;B(m,n) G [0,1], (12) 

m=l n=m+l 

Note that because D can take positive and negative 
values, might be in principle nullified not only by 
identical intensity traces, but also for those which are 
different but cross each other one or more times. By do¬ 
ing so we allow a pair of similar traces to be regarded as 
equal (i.e., tolerance in this case is better than if absolute 
value bars are introduced in Eq. (9)) but we risk that two 
very different traces are taken as equal (we therefore as¬ 


sume the probability for this to be small). This could be 
regarded as a drawback of this similarity measure and its 
significance depends in reality on the nature of the ana¬ 
lyzed data. However, the absolute value bars in Eq. (10) 
ensure that little non-zero values of will accumulate 
when performing the angular integrals (or sums), wors¬ 
ening the overall symmetry. To finalize, we stress here 
that we tried six different ways (not specified here) of 
obtaining symmetry indices and the one presented here 
proofed (unlike all others) to give systematically a good 
qualitative correspondence for all situations with differ¬ 
ent masks and input peak intensities. We recall that the 
purpose of is the one of obtaining a global feature of a 
beam profile with complex structure. 
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